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Turbulence modulation by inertial-range-size, neutrally-buoyant particles is investigated experi¬ 
mentally in a von Karrnan flow. Increasing the particle volume fraction <f> v , maintaining constant 
impellers Reynolds number attenuates the fluid turbulence. The inertial-range energy transfer rate 
decreases as a $ v , suggesting that only particles located on a surface affect the flow. Small-scale 
turbulent properties, such as structure functions or acceleration distribution, are unchanged. Fi¬ 
nally, measurements hint at the existence of a transition between two different regimes occurring 
when the average distance between large particles is of the order of the thickness of their boundary 
layers. 


Introducing impurities in a developed turbulent flow has drastic effects on the mechanisms of energy 
transfer and dissipation. A minute amount of polymer additives causes for instance drag reduction. ^ In 
such visco-elastic fluids, the coupling between the local flow and the polymers stretching can be modeled 
in order to quantify the exchanges between kinetic and elastic energies and to interpret turbulent drag 
reduction as a suppression of large velocity gradientsP The mechanisms at play in turbulent suspensions of 
finite-size spherical particles are much more intricate and turbulence can either be enhanced or suppressedP 
In wall flows, transition to turbulence can be either hindered or facilitated depending on the size and volume 
fraction of particlesP^ In the case of developed turbulence, depending on the size of impurities, on their 
density ratio, and volume fraction, both turbulent kinetic energy and dissipation rate can either increase 
or decrease! 6 ‘ There is no thorough understanding of turbulence modulation by finite-size particles. It 
is only in the asymptotics of very small and very dilute particles that one can satisfactorily model the 
energy budget between the fluid and the dispersed phase. Away from these limits, either finite-size effects 
preclude an explicit formulation of the forces exerted on the particles or large concentrations require the 
understanding of possible collective effects. 

Finite-size effects are usually accounted for in terms of Faxen corrections to the Stokes equation governing 
point-particle dynamics!^ However, it was observed 10 that this approximation works only for particle whose 
diameters is below « 4 t 7, where rj designates the turbulent Kolmogorov dissipative scale. The finiteness of the 
particle Reynolds number becomes important at larger sizes and explicit models would require solving the 
full non-linear Navier-Stokes equation with the proper boundary conditions at the particle surface. A first 
step in tackling finite-size effects consisted in characterizing the statistics of various dynamical quantities, 
such as velocity, acceleration, rotation rate, for large neutrally buoyant particles in developed turbulent 
flows P 1 13 One of the main findings is that average particle dynamical properties can be predicted from 
usual turbulent dimensional analysis at scales of the order of the particle diameter. More recently, much 
effort has been devoted to the characterization of the fluid flow in the vicinity of a single particleP^*® i n 
particular, it was found that it is surrounded by a boundary layer whose average thickness is of the order 
of its diameter. At larger distances, the influence of the particle onto the carrier fluid is negligible. Also, 
kinetic energy dissipation is st rongl y increased at the particle surface (in its viscous boundary sublayer) and 
weakly diminished in its wake. 15 11 

All these observations for neutrally spherical particles were obtained in the limit where they are very 
dilute, so that individual disturbances of the carrier never pile up nor interfere. At higher concentrations, 
large particles interact together through either collisions or perturbations of the carrier flow. This “four-way 
coupling” occurs when particles are at distances less than the thickness of their boundary layer. In principle, 
if the particles are uniformly distributed in the flow with a volume fraction $ v < 1, the effects of two-particles 
interactions are expected to be oc <f> 2 , and are thus negligible compared to single-particle contributions which 
are oc <F V . Several works indicate that such a naive picture cannot be true. First, neutrally-buoyant particles 
seem to distribute in a non uniform manner; they cluster near walls 18 19 or correlate with the large scales 
of the carrier flowP3 It is thus probable that local particle density has strong fluctuations. Second, large 
particles that approach very close to each other might spend some time together. This is suggested for 
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instance in the simulations of Ref. cm and could be due to strong dissipative mechanisms occurring when 
particles meet. This effect could lead to the creation of particle clumps, as for instance observed for heavy 
small-size elastic particles.^ Such collective phenomena, if exist, might be crucial in the understanding of 
turbulence modulation by finite-size particles. 

We present here experimental results on the influence of finite-size, neutrally buoyant particles on a 
turbulent von Karman flow. We find that the most noticeable effects induced by the dispersed phase is 
a continuous decrease of the turbulent kinetic energy and of the average inertial-range transfer rate when 
increasing the particle volume fraction while keeping constant the rotation speed of the impellers. We also 
obtain evidence that the small-scale turbulent properties of the fluid flow are unchanged by the presence 
of the particles. In particular, our results show surprisingly that in the presence of particles, second-order 
statistics of velocity increments match those of unladen flows. Finally, our results support the existence of 
a transition between two different regimes which occurs when the average distance between large particles 
is of the order of the thickness of their boundary layers. 

We consider a von Karman water flow maintained in a developed turbulent state (with R\ rs 300) by 
two impellers of diameter d = 28 cm that counter-rotate at an approximately constant torque. Specifically, 
we adjust the power of the motors at the beginning of each experiment in order to assert the same average 
rotational frequency of the impellers for all measurements. The enclosure has an octagonal cylindrical shape 
with dimensions are 40 x 38 x 38 cm. The flow field is analyzed by tracking the temporal evolution of 
fluorescent tracer particles (D rs 107 ^m ) using three cameras (Phantom V10, manufactured by Vision 
Research Inc., Wayne, USA) at 2900 fps. The measurement volume of approximate diameter 8cm in the 
center of the tank is illuminated by a 100W laser beam. We use the weighted averaging algorithm for two- 
dimensional tracer finding and a “three-frames minimum acceleration” for tracking their three-dimensional 
trajectories? 23 The finite-size particles are super-absorbent polymer spheres whose optical index and mass 
density match those of water. These particles have diameters of the order of D p rs 90?7 ~ L/ 9, where rj 
and L are the dissipative and integral length scales of the flow without particles. Two thin grids situated 
at 10 cm from the impellers prevent particles from colliding with them. This setup has been already used in 
Ref. \M and [2H an d more details can be found therein. We study here the fluid flow characteristics varying 
the volume fraction of large particles from <I> V = 0 to 10%. Each time the particle load is changed, the power 
of the impellers is adjusted to maintain constant their rotation frequency at /i mp = 0.79 Hz (±1.5%), and 
thus the associated Reynolds number at Re = d 2 /i mp /^ ~ 62 000. The different setups are summarized in 
Table H 


4>v 

0 

0 .1% 

0.3% 

0 .8% 

2 % 

4% 
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/imp (Hz) 

0.79 

0.79 

0.79 

0.79 

0.80 

0.80 

0.80 

w rms (ms -1 ) 

0.0796 

0.0797 

0.0771 

0.0752 

0.0766 

0.0578 

0.0466 

e (m 2 s -3 ) 

0.0058 

0.0055 

0.0053 

0.0048 

0.0051 

0.0032 

0.0012 

V (Mm) 

115 

117 

118 

121 

119 

133 

183 

r v (ms) 

13 

14 

14 

15 

14 

18 

33 

R\ 

321 

335 

319 

319 

316 

228 

143 


TABLE I. Characteristics of the various experiments, 4> v volume fraction of the large particles; /i mp average 
frequency of the impellers; M rma root-mean squared velocity (averaged over components); e average kinetic energy 
transfer rate (measured from second-order structure function; see text); r] = i' 3 ' /4 /g 1 / 4 Kolmogorov dissipative scale; 
t v = v 1 / 2 /e 1 / 2 associated turnover time; R\ Taylor-microscale Reynolds number. 

We start by considering the effect of the large particles on the intensity of turbulent fluctuations. The 
Eulerian mean velocity profile ( u{x )) is obtained for the different values of <J> V by binning positions and 
averaging velocities over all movies. We then define velocity fluctuations as instantaneous deviations of the 
tracers velocity from this average field. Figure [lja) represents the variance u 2 ms of the axial (along the axis 
of symmetry) and transverse components of these turbulent velocity fluctuations as a function of the volume 
fraction of large particles. One clearly observes that the turbulent kinetic energy decreases when the particle 
load increases, as already observed in a different turbulent flow. ' Data show anisotropy, which persists for 
all $ v : the root-mean square values of axial velocity components are around 60-70% of the transverse ones. 
Measurements are in good agreement with w 2 ms ($ v ) rs w 2 ms (0)(l — 3.75 <hv' /3 ). For transverse components, 
one point, corresponding to <J> V rs 2 % clearly deviates from this form. As we will see later, this could be due 
to difficulties in estimating the mean flow in that case. The fitting form predicts that it rms = 0 for <f> v > 14%. 
This could correspond to a critical volume fraction above which turbulence is completely extinguished by 
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the particles. 



FIG. 1. (a) Variance of the turbulent velocity fluctuations u 2 ms = ((it* — (u;}) 2 ) as a function of the particle 
volume fraction; The dashed lines are «r ms (^v) = u?ms(0)(l — 3.75 <f > i/ 3 ). (b) Discrepancy of the average inertial- 
range kinetic energy transfer rate Ae = 1 — e($ v )/e(0) as a function of the particle volume fraction <f> v ; the dashed 
line is Ae = 3.75 The rate e is estimated from the inertial-range average value of the compensated second- 

order longitudinal velocity structure functions S 2 (r)/(C 2 r 2 ^ 3 ) using for the Kolmogorov constant C 2 = 2.1. Inset: 
Dissipation rate obtained from the second-order Eulerian longitudinal structure functions for the various volume 
fractions (the upper curve corresponds to the case with no particles and <f> v increases from top to bottom). 


The inertial-range kinetic energy transfer rate e follows a similar trend. Figure mb) shows that the 

2/3 —^ 

discrepancy Ae = 1 — e(<fr v )/e(0) is proportional to dv , like for the turbulent kinetic energy. Even at 
small values of the volume fraction, the attenuation of turbulence does not show a linear behavior oc <f> v 
(proportional to the total number of particles) that is expected if the effect of particles consists in the 
superposition of non-interacting individual perturbations. An attenuation oc <f>v^ 3 suggests that collective 
effects and interactions between particles are at play. Such a scaling indicates that only a fraction of the 
particles have an important impact on the flow. The power two-third can be interpreted as if these active 
particles are located on a surface rather than in the full volume of the experiment. Several heuristic scenario 
can lead to such a law. First, we can imagine that particles form a given number of aggregates where all 
fluid flow perturbations occur in the boundary layers of the outer edge. We did not detect such aggregates 
in the observation volume but we cannot exclude that they are located in regions of the flow that we do not 
monitor (for instance close to the grids near the impellers). Recent work found that particles are likely to 
cluster far from the center of the von Karman flowP^ A more likely scenario relies on a kind of shielding 
effect where the kinetic energy is prevented from reaching the center of the flow by the particles situated 
at the periphery of the experiment. A full discrimination between these possibilities requires systematically 
scanning the entire experiment in order to measure the inhomogeneities in both the particle distribution 
and the fluid turbulence characteristics. 

We next move to small-scale turbulent properties. A striking result is that they seem unaffected by 
the presence of the particles once scaled with the corresponding transfer rate e, large-scale kinetic energy 
content u' 2 ms . and associated Reynolds number R\ = u 2 m s y/lh/{ve). Figure [i|a) shows the second-order 
longitudinal Eulerian structure function of the turbulent velocity fluctuations for the various particle loading 
considered. All curves collapse in the inertial and sub-inertial ranges after they are represented in dissipative- 
scale units. In particular there are no noticeable effects (up to noise) at separations of the order of the 
particles diameter for <f> v up to 10%. These results seem to contradict the kinetic energy spectra obtained 
in numerical simulations 21 25 that display an enhancement of the energy contained at scales smaller than 
the particle size and a depletion at larger scales. These spectra are obtained from a Fourier transform of the 
full domain encompassing both the fluid and the particles. There are still questions whether such effects on 
the spectrum are inherent to the perturbed fluid flow or come from the particle velocity field. In our case, 
the use of Lagrangian tracers to perform Eulerian measurements isolate the fluid statistics from such pitfalls 
and suggest that there is no intrinsic modifications of the distribution of turbulent velocity increments by 
the large particles. 

Statistics along tracer trajectories seem also rather insensitive to the presence of large particles. Fig- 
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FIG. 2. (a) Second-order Eulerian longitudinal structure functions S 2 (r) in dissipative-range units obtained for 

the various values of the volume fraction 4> v ; The solid line corresponds to a behavior oc r 2//3 ; The two vertical 
dashed lines shows the range 54 < Dp/rj < 88 obtained from the maximal and minimal values of r/. (b) Second-order 
Lagrangian structure function S^ s (t) in dissipative-scale units. Inset: logarithmic derivative (dlog log S^ s ) 

of the fourth-order Lagrangian structure function with respect to the second-order. 


ure[2|b) shows measurements of the Lagrangian structure functions Sp ag = ([ Ui(t + r) — Ui(t)] p ), where Ui 
are components of the tracer velocity and the average is over the ensemble of trajectories. One observes 
that for p = 2, the collapse is not as good as for Eulerian statistics. However, Lagrangian structure func¬ 
tions are known to display very intermittent properties and the observed discrepancies are comparable to 
those obtained when comparing unladen flows at different Reynolds numbersP^We see from Tab.jljthat our 
Reynolds number decreases when the particle volume fraction increases: R\ varies by more than a factor 
two when goes from 0 to 10%. This lack of universality is not present when considering only the scaling 
properties through the logarithmic derivatives of the Lagrangian structure functions. For that reason, we 
follow Ref.[26]by showing in the inset of Fig.^b) the quantity (dlog ^^/(dlogS'Jf 8 ). Results from different 
volume fractions collapse to the same universal curve. 

The statistics of the fluid acceleration seem also not affected by the large particles. Figure |3](a) shows the 
density functions of the acceleration components after they are normalized by their respective second- 
order moments. Up to possible slight variations in the tails that could correspond to the change in Reynolds 
number at varying the fraction of large particles, the different curves seem to collapse reasonably well. It is 
known that the moments of these distributions do not reflect the actual value of the acceleration variance 
because of data filtering. To measure (a 2 ) we follow the method of Voth et consisting in extrapolating 
the variance when decreasing the filter size to zero. The resulting values of ao = (a 2 )zA/ 2 /e 3 / 2 are shown 
in the inset of Fig. §b) as a function of the Reynolds number of the flow. One clearly observes that 
the decrease of ao when <f> v increases relates to the change in Reynolds number. Finally, the acceleration 
autocorrelations represented in the main panel of Fig. |3|b) are not much affected by the particle presence. 
The slight variations are again due to the change in Reynolds number induced by the particles. 

All the quantities considered above are either reduced as $v^ 3 or not affected by the large particles. They 
relate to the turbulent fluctuations from which we have removed the mean flow. We actually find that the 
mean flow itself is strongly affected by the large-particle load. To quantify its variations as a function of the 
particle volume fraction, we have computed its kinetic energy large-particle load. To quantify its variations 
as a function of the particle volume fraction, we have computed its kinetic energy F(<f> v ) = f |(it(:r))| 2 d 3 :r. 
Figure |4](a) shows that F($ v ) remains almost constant up to ss 3% where it suddenly drops down. Above 
this critical value, the mean flow is in a different state that contains approximately twice less energy. The 
transition occurs roughly when the typical distance between particles becomes less than their diameter. All 
particles are thus necessarily feeling the presence of their neighbors through the flow modifications. This 
apparently has a drastic impact on the average motions as can be seen when comparing Figs. [4](b) and (c). 
They represents the variations of the energy and of the (x, z) components of the mean flow in the same 
slice of the measurement volume for <f> v = 0.8% (b) and 4% (c). Not only is the energy content strongly 
modified but also the fine structures and the directions of the mean flow. Panel (b) is representative of 
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FIG. 3. (a) Normalized probability density functions of acceleration components for various values of the volume 

fraction 4> v as labeled; The black dashed line is the approximative form obtained when assuming that the acceleration 
amplitude is log-normal. 2 ^ (b) Temporal autocorrelation of the acceleration components for the same values of $ v - 
Inset: Normalized variance of the acceleration as a function of the Reynolds number R\; Each symbol is a different 
experiment and the dashed line is the fit proposed by Hill 29 for the intermittent dependence of the acceleration 
variance as a function of R\. 
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FIG. 4. (a) Mean flow kinetic energy -E(<f> v ) contained in the measurement volume as a function of the large-particle 

volume function, (b) and (c) Slice of the mean flow in a central part of the measurement volume for <f> v = 0.8% and 
4%, respectively (the x direction is along the impellers axis of rotation). The colored background shows the velocity 
squared modulus (normalized in each case to its mean; increasing values from white to black) and the arrows shows 
the ( x, z) components of the velocity field. 


all volume fractions below 0.8% while panel (c) is also representative of <!> v = 10%. For <h v = 2%, the 
mean flow obtained when time-averaging is very different from these two states. We however suspect it to 
be non stationary because several movies show extremely strong deviations from the mean. Unfortunately, 
our experiment was conceived to focus on small-scale turbulent properties and not designed to accurately 
measure large scale variations. A better understanding of what is happening for this volume fraction would 
require revisiting our setup in order to focus on the large scales. This would also allow us a more accurate 
study of a possible phase transition for $ v of the order of a few percents. 

To conclude, let us stress again our main findings on the modification of a turbulent von Karman flow by 
finite-size neutrally-buoyant spherical particles. The most surprising result is that inertial-range and small- 
scale turbulent features are unchanged in the presence of particles for volume fractions up to d> v = 10%. 
There is no signature of the particle size and its associated timescale neither on the scaling properties of the 
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fluid turbulent velocity fluctuations nor on the statistics of the acceleration. The only noticeable effects are 
on global quantities. In our specific setup and when the impeller rotation rate is maintained, the turbulent 
kinetic energy and the inertial-range transfer rate of the bulk flow are continuously decreasing when the 
particle volume fraction is increased. However, our setup does not monitor the energy injection rate. We 
are thus unable to discriminate between a decrease of the power needed to maintain the impellers rotation 
speed (i.e. a sort of drag reduction) and a redistribution of the energy dissipation to the regions closer to 
the impellers (outside of our observation volume). A better understanding would require a more accurate 
handling of the impeller torque and thorough measurements of the fluid velocity in the full experimental 
domain. Such measures could clarify the possible phase transition at 4> v ~ 3%. In particular, they would 
allow us to draw differences between what is due to the specificities of von Karman inhomogeneous flow 
and a possible universal change at such volume fractions in the coupling between the turbulent flow and the 
particles. 
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